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ABSTRACT 


A  computer  program  has  been  written  to  calculate 
the  heat  transfer  on  the  windward  side  of  a  delta  wing 
in  hypersonic  flow.  The  program  utilizes  the  small 
crossflow  theory  of  Vaglio-Laurin.  Heat  transfer  values 
are  calculated  for  both  the  laminar  and  turbulent  boundary 
layer  cases,  and  results  are  compared  with  experimental 
data  on  delta  wings  available  in  the  literature. 

ADMINISTRATIVE  INFORMATION 

This  work  was  authorized  and  funded  jointly  by  the  Naval  Air  Systems 
Command  and  National  Aeronautics  and  Space  Adminstration  Marshall  Space 
Flight  Center  under  Naval  Air  Systems  Command  Task  A370538B  and  NASA  Purchase 
Request  H58564A. 


INTRODUCTION 

The  objective  of  this  work  was  to  develop  a  computer  program  to 
predict  the  heat  transfer  characteristics  on  the  windward  side  of  a  delta 
wing  at  hypersonic  speeds.  The  program  could  then  be  used  to  determine  the 
thermal  environment  of  the  configuration. 

1  2* 

The  method  of  approach  was  based  on  the  work  of  Vaglio-Laurin,  * 
who  demonstrated  that  on  three-dimensional  bodies  with  highly  cooled  surfaces 
and  moderate  Mach  numbers  at  the  outer  edge  of  the  boundary  layer,  crossflow 
terms  in  the  boundary  layer  were  negligible--even  for  large  transverse 
pressure  gradients.  With  this  assumption,  heat  transfer  characteristics  on 
the  three-dimensional  body  could  be  found  using  formulas  derived  for  ax- 
isymmetric  (no  crossflow)  bodies  provided  that  the  pressure  distribution  on 
the  three-dimensional  body  be  known. 


^Vaglio-Laurin,  R. ,  "Laminar  Heat  Transfer  on  Three-Dimensional  Blunt 
Nosed  Bodies  in  Hypersonic  Flow,"  America  Rocket  Society  Journal  New  York, 
p.  123-129  (Feb  1959). 

2 

Vaglio-Laurin,  R.,  "Turbulent  Heat  Transfer  on  Blunt-Nosed  Bodies  in 
Two-Dimensional  and  General  Three-Dimensional  Hypersonic  Flow,"  Journal  Aero/ 
Space  Science  New  York,  p.  27-36  (Jan  1960). 

* 

A  complete  list  of  References  is  listed  on  page  35. 
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The  calculation  is  in  two  parts.  First,  the  inviscid  streamline 
pattern  and  thermodynamic  property  distribution  in  streamline  coordinates 
are  found  on  the  body  surface.  Second,  at  each  calculated  point  on  the  body 
surface,  analogy  is  made  with  various  axisymmetric  heat  transfer  equations, 
and  the  heat  transfer  characteristics  are  determined.  The  method  of  Lees 

and  Vaglio-Laurin'*'  are  used  for  laminar  flow;  a  streamline  divergence  tur- 

4  5 

bulence  flow  method  and  a  modified  method  by  Rose,  Probstein,  and  Adams 

are  used  for  turbulent  flow. 

The  actual  calculation  begins  from  the  stagnation  point  on  the 
spherical  nose,  proceeds  to  the  cylindrical  leading  edge,  and  goes  on  to  the 
flat  plate  regions.  Appropriate  integrals  of  thermodynamic  properties  are 
calculated  along  the  streamlines  to  permit  heat  transfer  analysis.  Various 
pressure  distributions  are  used  in  each  region.  The  boundary  layer  entropy 
values  used  herein  were  determined  from  the  normal  shock  at  the  nose. 

The  report  presents  an  analytic  development  of  the  small  crossflow 
approximation  as  derived  by  Vaglio-Laurin,  gives  the  specific  governing 
equations  for  the  streamline  pattern  on  a  delta  wing,  and  presents  the 
pressure  distribution  expressions  that  apply  to  this  geometry.  Then  the 
thermodynamic  variables  required  in  the  analysis  are  discussed  along  with 
the  governing  heat  transfer  equations  that  use  these  variables.  Results 
obtained  by  using  this  theory  are  compared  with  experimental  results  avail¬ 
able  in  the  literature  for  blunt-nosed  hypersonic  layers.  These  comparsions 
include  calculated  and  experimentally  determined  pressures,  resulting  stream¬ 
line  patterns,  and  computed  and  measured  heat  transfer  values. 


3 

Lees,  L.,  "Laminar  Heat  Transfer  over  Blunt-Nosed  Bodies  at  Hypersonic 
Flight  Speeds,"  America  Rocket  Society  Journal  New  York,  p.  259-269  (Apr  1956). 

4 

Ratliff,  A.W.  et  al.,  "Analysis  of  Heating  Rates  and  Forces  on  Bodies 
Subject  to  Rocket  Exhaust  Plume  Impingement,"  Huntsville,  Ala.,  Lockheed 
Missiles  §  Space  Co.,  LM  SC/HREC  A79  1230,  HREC  1150-1,  Contract  NAS8-21150 
(Mar  1968). 

^Tai,  T.C.,  "Laminar  and  Turbulent  Convective  Heat  Transfer  over  Bodies 
at  an  Angle  of  Attack."  Virginia  Polytechnic  Institute,  Thesis,  Blacksburg, 

Va.  (Oct  1968). 


2 


ANALYTICAL  DEVELOPMENT 


The  equations  for  the  steady, laminar  boundary  layer  flow  of  a  homo¬ 
geneous  gas  in  the  absence  of  body  forces  and  heat  sources  may  be  written 
in  an  orthogonal  (x^,  x^,  x^)  coordinate  system  as  (in  this  section  only  hf 
refers  to  local  enthalpy  and  H  to  stagnation  enthalpy) 


2 

These  results  are  taken  directly  from  Vaglio-Laurin.  Choosing  x^ 
coincident  with  the  inviscid-f low  streamline  at  the  outer  edge  of  the 
boundary  layer  and  x^  normal  to  the  surface,  one  has  as  boundary  conditions 

x3  =  0  V1  =  V2  =  V3  =  0  ;  H  =  Hw  (2a) 


x 


3 


OO 


V2  -  0  ;  H  =  He 


(2b) 
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For  the  case  of  highly  cooled  walls  and  V^e/(Hg  -  Hw)  <<  1 


-*•  0 


Vle  p 


(3) 


As  discussed  in  Reference  1,  it  can  also  be  shown  that  under  the  subject 
conditions,  V2  =  0  (inside  of  the  boundary  layer).  Noting  that  h^  =  1, 
the  governing  equations  become 


with  the  boundary  conditions 


x_  =  0  V.  =  V_  =  0  ;  H  =  H  (5a) 

3  13’  w  v  ' 

x3  V1  =  Vle=  Ue  H  =  He  ™ 

Examination  of  the  preceding  equations  shows  that  they  become 
formally  identical  to  those  governing  the  boundary  layer  over  an  axisymmetric 
body  at  zero  angle  of  attack,  provided  that  be  replaced  by  the  radial 
coordinate  r  and  distance  be  measured  along  respective  streamlines.  The 
same  equations  that  govern  the  heat  transfer  on  an  axisymmetric  body  at 
zero  angle  of  attack  may  also  be  used  in  the  present  nonaxisymmetric  case, 
provided  that  r  be  replaced  by  in  the  axisymmetric  heat  transfer  equations. 

Three  regions  on  the  delta  wing  are  under  consideration:  the  spherical 
nose  cap,  the  cylindrical  leading  edges,  and  the  flat  plate  lower  area. 
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Under  the  assumptions  made  in  this  study,  all  boundary  layer  streamlines 
go  from  the  geometric  stagnation  point  on  the  spherical  nose,  to  the 
cylindrical  leading  edge,  to  the  flat  plate  region.  Thus  it  is  necessary 
to  have  the  governing  equations  in  a  form  appropriate  for  calculating  the 
streamline  motion  and  the  corresponding  values  of  h^  in  each  region. 

GOVERNING  DELTA  WING  STREAMLINE  EQUATIONS 

Although  the  directions  of  the  streamlines  on  the  spherical  nose 
are  known  in  advance  to  be  great  circles  through  the  stagnation  point,  ad¬ 
ditional  variables  must  still  be  calculated  in  order  to  determine  h^  at 
each  location.  The  equations  used  directly  to  do  this  are  available  in 
Reference  5  and  are  not  repeated  here. 

The  following  equations  for  the  cylindrical  and  flat  plate  areas 
were  derived  by  Tai0  as  an  extension  of  the  work  reported  in  Reference  5. 
The  cylindrical  and  flat  plate  regions  are  shown  in  Figures  1  and  2. 

The  governing  equations  for  the  cylindrical  region  are 


Dx. 

— -  =  cos  6 


D<}>b  sin  0b 

DS  R 
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D0, 
_ b 

DS 


9P_ 

3xb 


COS0b  9P  \ 
Ro  **b/ 


^2  =  i_^b 

DS  9  & 


(6a) 


(6b) 


(6c) 


g 

Tai,  T.C.,  "Streamline  Geometry  and  Equivalent  Radius  over  a  Flat 
Delta  Wing  with  Cyclindrical  Leading  Edge  at  Angles  of  Attack,"  NSRDC  Report 
3675  (Oct  1971). 
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The  line  <j)  equals  zero  corresponds  to  the  zero  angle  of  attack  cylindrical 
stagnation  line. 

The  flat  underside  of  the  wing  is  governed  by  the  following  equations 


Dxc 

=  cos6c 


(7a) 


D<f)  sin0 


DS 


(7b) 


D0 
_ c 

DS 


yM2P 
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9P 
c  3x 


COS0 


c  9P 
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sin0 


(7c) 


Dh, 


DS  9  e 


30  h„cos0 

c  2  c 


(7  d) 
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It  is  also  important  to  know  the  joining  conditions  between  the  sphere 
and  the  cylinder  and  between  the  cylinder  and  the  flat  plate  regions. 

Figure  3  is  a  diagram  of  the  sphere  region.  The  streamline  that  is  being 
analyzed  goes  off  the  sphere  and  onto  the  cylinder  when 

S/Ro  >  f  (8) 

where  S  is  the  streamline  length 
Rq  is  the  sphere  radius  and 

f  (see  Figure  3)  is  as  defined  in  Equation  (10b). 

The  initial  conditions  on  the  cylinder  are  determined  from  the  final 
conditions  on  the  sphere  and  from  the  coordinates  and  equations  of  the 
cylindrical  region. 

=  0  (9a) 

l 

V  =  i  -  b  (9b) 

l 
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(9c) 


(9d) 


E  and  b  are  spherical  dimensions  shown  in  Figure  3.  The  unknowns  in  the 
preceding  development,  namely,  E  and  f,  may  be  readily  determined  from  the 
law  of  sines  and  cosines  as 


when  <f>b 


E  =  cos 


(-cos  F  cos  B  +  sin  F  sin  B  cos  e) 


F  =  sin 


-1 


(' 


sin  f 


sm  e 
sin  E 


) 


(10a) 

(10b) 


The  streamline  goes  off  the  cylinder  and  onto  the  flat  plate  region 
>  90  degrees.  The  initial  conditions  on  the  flat  plate  are 


c. 

l 


(11a) 


d>  =0 
c . 


(lib) 


l 


(11c) 


(Hd) 


DELTA  WING  PRESSURE  DISTRIBUTION  EQUATIONS 

The  pressure  distribution  must  be  known  on  each  region  of  the  delta 
wing  in  order  to  apply  the  governing  streamline  equations.  This  section 
presents  various  expressions  for  pressure  distribution. 
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In  the  spherical  nose  region,  a  modified  Newtonian  pressure  distri¬ 
bution  is  used  in  the  numerical  calculations 


2 

cos 


(12) 


This  pressure  distribution  is  a  function  only  of  distance  from  the  stagnation 

point.  This  formula  can  be  expected  to  be  in  greatest  error  where  the  sphere 

merges  to  the  cylinder  when  the  delta  wing  is  at  low  angles  of  attack. 

For  the  cylindrical  leading  edges  and  flat  plate  area  of  the  delta 

wing,  the  pressure  equations  are  found  directly  from  the  theory  presented  in 

7 

Buck  and  McLaughlin.  The  pressure  equation  for  the  cylindrical  region  is 


y-  =  C  (  0.320  +  0.455  coscJ>£  +  0.195  cos2<|>£ 

roo  ' 

+  0.035  cos3(j>^  -  0.005  cos4<j^\  (13a) 


C 


1 


(Y+l)Mro2  cos2Ae 

_x_ 

Y-l 

Y+l 

2 

_2yM2cos2Ae  -(y-l)_ 

(13b) 


where  is  the  effective  sweep  angle  of  the  wing  as  defined  in  Reference  7, 
and  <f>k  is  a  cylindrical  angle  measured  with  respect  to  the  effective  leading- 
edge  stagnation  line. 

Note  that  the  previous  expression  for  pressure  distribution  is  a 
coordinate  function  solely  of  <J>^.  Equations  (13)  are  based  on  correlated 
cylinder  data.  The  cylindrical  pressure  equations  include  corrections  for 
the  angle  of  attack  and  sweep,  the  neglect  of  the  nose  region,  and  the 
neglect  of  the  afterbody  effect  of  the  flat  plate  area,  that  is,  the  cylin¬ 
der  does  not  continue  through  after  the  point  where  =  90  degrees.  This 


7 

Buck,  M.L.  and  McLaughlin,  E.J.,  "A  Technique  for  Predicting  Pressure 
Distributions  and  Aerodynamic  Force  Coefficients  for  Hypersonic  Winged  Re- 
Entry  Vehicles  (U),"  Wright-Patterson  Air  Force  Base,  Aeronautical  Systems 
Division,  TDR-63-522  (Aug  1963).  CONFIDENTIAL 
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last  mentioned  factor  is  considered  to  account  for  the  small  pressure  jump 
between  the  preceding  cylindrical  formula  and  the  flat  plate  formula  shown 
as  follows. 

The  expression  for  pressure  on  the  windward  flat  plate  region  is 


(■ 


+  BC 

Y 


2/3 

2  2 
M  cos^A 

CO 


(14) 


The  first  right-hand-side  term  is  a  Newtonian-type  expression.  K  is  a 
function  of  angle  of  attack;  K  versus  a  is  given  in  Reference  7  where  the¬ 
oretical  results  are  compared  with  flat  plate  experimental  data.  The  second 
small  additional  term  is  what  creates  the  pressure  variation  on  the  flat 
plate  region.  This  term  is  based  on  an  analogous  blast-wave-type  effect 
across  the  blunt  cylindrical  leading  edge.  The  most  important  fact  to  note 
in  the  second  expression  is  that  the  pressure  change  produced  is  normal  to 
the  leading  edge,  x^-direction.  The  constants  appearing  in  this  expression 
are  given  in  Reference  7.  The  present  analysis  does  not  consider  boundary 
layer-displacement  effects  on  the  pressure. 

It  is  difficult  to  determine  the  value  of  A,  the  shock  standoff 
distance  along  the  leading  edge.  In  the  computer  program,  it  was  set  equal 
to  the  standoff  distance  in  front  of  the  spherical  nose.  Pressure  values 
determined  by  other  theories  may  be  readily  substituted  into  the  program 
in  place  of  the  preceding  results,  provided  that  the  derivatives  of  pressure 
can  be  explicitly  evaluated.  Finally,  it  should  be  noted  that  in  the  vi¬ 
cinity  of  the  centerline  when 


tt/2  -  A  -  <{>c 
tt/2  -  A 


<  0.06 


(15) 


the  pressure  derivatives  have  been  set  equal  to  zero  in  the  program. 

THERMODYNAMIC  VARIABLES 

The  upstream  flow  values  Mra,  Tq,  Pq,  y  ,  Pr,  Cp,  and  y  are  known. 

The  boundary  layer  is  assumed  to  be  formed  of  upstream  flow  that  has  passed 
through  a  normal  shock  wave  in  front  of  the  spherical  nose  cap;  this  point 
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will  be  discussed  further  in  a  later  section.  The  stagnation  pressure 
of  the  flow  entering  the  boundary  layer  is  therefore 


Y  1 


(Y+l)Mco2 

Y-l 

Y+l 

_ (Y+IJM^2  -2  _ 

_2yMco2  -  (y-l)_ 

(16) 


The  value  of  pressure  in  the  boundary  layer  is  taken  as  equal  to  the 
inviscid  pressure  at  the  edge  of  the  boundary  layer.  The  inviscid  pressure 
at  each  computed  point  is  determined  from  Newtonian  or  other  theories  that 
will  be  described.  The  local  Mach  number  on  the  surface  of  the  body  (for 
flow  that  has  passed  through  the  shock)  is  then  found  from  the  isentropic 
flow  equation 


(17) 


To  determine  the  heat  transfer,  various  thermodynamic  values  must 
first  be  determined  at  infinity,  at  the  edge  of  the  boundary  layer,  at  re¬ 
ference  conditions,  and  at  the  wall.  The  four  heat  transfer  equations  used 
in  the  analysis  can  be  written  in  terms  of  these  values. 

The  flow  values  in  each  region  are  readily  calculated.  At  infinity 
upstream,  the  necessary  flow  values  are 


P 

oo 


T 

p 

oo 

P«>  =  RrT 
G  00 


(18a) 


(18b) 


(18c) 


11 


English  engineering  units  are  used  with  time  in  hours,  and  therefore  the 
calculations  of  velocity  values  are  performed  in  the  program  with  a  con¬ 
version  factor  of  3600.  Finally,  the  reference  conditions  are 

8 

Shapiro,  A.H.,  "The  Dynamics  and  Thermodynamics  of  Compressible  Fluid 
Flow."  Vol.  1.  Ronald  Press,  New  York,  (1953). 

9 

National  Advisory  Committee  for  Aeronautics,  "Equations,  Tables,  and 
Charts  for  Compressible  Flow,"  NACA  Report  1135  (1953). 

^Schlichting,  H.,  "Boundary-Layer  Theory,"  Sixth  Edition  McGraw-Hill, 
Inc.,  New  York  (1968). 
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(20a) 


T  =  0 . 5T  +0.5T  +0.22/T  -  T  \ 

r  e  w  (  aw  e  / 


where 


T  +  110\ 
-= - 

T  +  110/ 
e  ' 


P 


r 


P 

e 


Vr 


T 

aw 


T  + 
e 


(rc  »2  ^  RGTe) 
2Cp  *  778 


(20b) 


(20c) 


(20d) 


(20e) 


1/2  1/3 

and  rc  is  the  recovery  factor  set  equal  to  in  laminar  flow  and  to  P^, 

in  turbulent  flow. 

The  wall  temperature  T  is  read  into  the  program  at  a  constant  value. 
11  w 

The  Eckert  value  of  reference  temperature  is  used  in  the  majority  of  the 
calculations  for  obtaining  the  reference  values.  It  is  felt  that  best 
overall  results  will  be  obtained  by  using  a  compromise  temperature  value  as 
previously  described.  Vaglio-Laurin,  however,  suggests  other  values  for 
reference  temperature,  and  some  calculations  are  subsequently  discussed 
using  this  approach. 

Enthalpy  terms  required  in  the  heat  transfer  expressions  are:  stag¬ 
nation  enthalpy  (Hq  =  .  Tq),  wall  enthalpy  (H^  =  .  Tw) ,  and  adiabatic 

wall  enthalpy  (H  =  C  .  T  ).  Heat  transfer  values  are  found  in  the 
K  aw  p  aw' 

formulas  to  be  presented  using  both  the  terms  (H  -H  )  and  (H  -  H  ) .  The 

0  w  aw  w 

^Rohsenow,  W.M.  and  H.Y.  Choi.,  "Heat,  Mass,  and  Momentum  Transfer," 
Prentice  Hall,  Englewood  Cliffs,  N.J.  (1961). 
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expression  (H  -  H  )  is  based  on  local  enthalpy  values  and  is  felt  to 
give  the  best  heat  transfer  results.  Finally,  it  should  be  noted  that  y  is 
common  between  all  subroutines  and  that  its  value  may  be  changed  in  the 
course  of  a  computation  to  represent  real  gas  effects.  (This  was  not  done 
in  any  of  the  present  calculations) 

HEAT  TRANSFER  EQUATIONS 

Utilizing  the  previously  found  thermodynamic  values,  the  equation 

3 

for  heat  transfer  using  the  Lees  method  for  laminar  flow  in  the  vicinity 
of  the  stagnation  point  is: 


wherd 
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The  preceding  heat  transfer  equations  are  derived  by  analogy  with 
usual  axisymmetric  formulas;  in  the  axisymmetric  formulas,  the  local  radius 
r  is  replaced  by  h^.  The  integrals  appearing  in  this  and  other  heat  trans¬ 
fer  methods  are  calculated  in  the  program  by  a  Simpson  rule  subroutine. 

The  integral  appearing  in  the  denominator  of  the  Lees  method  is  evaluated 
in  two  ways:  (1)  by  forcing  it  to  have  a  value  such  that  q^  =  qwQ  at  the 
stagnation  point,  and  (2)  by  initially  setting  it  equal  to  zero.  After  the 
first  few  integration  steps  (where  the  integration  method  is  still  unstable 
in  any  case),  both  approaches  produce  approximately  equivalent  results. 

In  the  Lees  method,  qwQ  is  the  stagnation  point  heat  transfer.  The 
value  of  this  quantity  is  calculated  from  a  Stanton  number  read  into  the 
program.  The  Lees  method  is  a  stagnation  region  solution  that  results  in 
a  decay  in  heat  transfer  from  the  stagnation  value  previously  described. 

The  method  therefore  has  greatest  validity  in  the  region  in  the  vicinity  of 
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stagnation  point  while  the  flow  is  still  laminar.  Other  techniques  are 
more  appropriate  for  points  further  downstream. 

A  second  technique  for  laminar  flow  is  given  in  the  Vaglio-Laurin 
paper  on  laminar  streamline  divergence.^ 
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The  preceding  expressions  are  again  based  on  the  -  r  analogy.  The 
calculation  in  the  program,  using  the  Vaglio-Laurin  laminar  method,  is  per¬ 
formed  using  two  values  for  enthalpy,  H  =  Haw  and  H  =  H  .  The  Vaglio-Laurin 
laminar  flow  method  is  to  be  used  in  laminar  flow  beginning  approximately 
20  degrees  from  the  stagnation  point. 

Two  techniques  are  included  in  the  program  for  the  case  of  turbulent 

4 

flow.  The  first  is  a  streamline  divergence  technique  where 
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It  is  also  derived  by  analogy  with  the  usual  heat  transfer  method  on  an 
axisymmetric  body  and  evaluated  here  using  H  =  and  H  =  H&w.  The  second 
technique  is  the  method  of  Rose  et  al .  as  modified  by  Tai^  to  include  pres¬ 
sure  gradients  effects.  This  formula  is 
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In  the  program,  results  are  determined  as  the  ratio  of  the  local  Stanton 
number  to  the  Stanton  number  at  the  stagnation  point. 

The  preceding  equations  along  with  expressions  for  pressure  distri¬ 
bution  presented  previously  are  used  in  the  computer  program  to  determine 
heat  transfer  results  along  calculated  streamline  paths.  The  computer 
program  is  not  included  in  this  report  (but  a  listing  will  be  kept  on  file). 
The  program  has  been  annotated  so  as  to  describe  the  operations  in  each 
section.  A  separate  sheet  is  maintained  to  describe  input  requirements 
and  control  statements.  Subroutine  form  has  been  used  to  perform  the 
calculations  in  the  sphere,  cylinder,  and  flat  plate  regions;  thermodynamic 
values  are  evaluated  in  groups  (boundary  layer  edge  values,  reference  values, 
etc.)  and  the  heat  transfer  equations  are  written  in  terms  of  this  common 
set  of  expressions.  An  attempt  has  been  made  to  make  the  program  easy  to 
follow;  comment  cards  have  been  included  where  appropriate,  and  variable 
names  have  been  made  to  correspond  to  the  spelling  of  physical  terms. 

RESULTS  AND  DISCUSSION 

Figures  4-11  compare  the  results  of  the  calculations,  using  the 
previously  described  theory  with  experimental  delta  wing  data  reported  in 
the  literature.  The  pressure  distributions  developed  in  the  previous 
section  are  first  compared  with  the  measured  pressure  values.  Then  after 
a  discussion  of  resulting  streamline  patterns,  comparsion  is  made  between 
computed  and  measured  heat  transfer  values. 

For  laminar  flow,  results  obtained  by  using  the  computer  program 
are  compared  with  experimental  National  Aeronautics  and  Space  Administration 
data  on  a  70-degree  slab,  delta  wing  (with  the  same  basic  geometry  as  used 
in  the  present  analysis).. in  hypersonic  flow  at  several  Mach  numbers  and 
various  angles  of  attack. 
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As  seen  in  Figure  4,  agreement  was  satisfactorily  obtained  for  the 

spherical  nose  region,  pressure  distribution.  For  the  most  part,  agreement 

12 

of  computed  and  measured  pressure  values  was  good  at  various  positions 
normal  to  the  leading  edge;  see  Figure  5.  At  high  angles  of  attack,  the 
flat  plate  part  of  the  delta  wing  will  strongly  influence  the  pressure  dis¬ 
tribution  on  the  cylindrical  leading  edge,  and  this  tends  to  make  the 
correlated  cylindrical  leading-edge  pressure  formula  in  Equation  (13a)  in¬ 
valid.  The  comparsion  of  experimental  and  theoretical  pressures  along  the 
centerline  of  the  wing  at  various  angles  of  attack  (Figure  6)  is  satisfactory. 

Figure  7  shows  computer-generated  streamline  patterns  on  the  same 
wing,  and  Figure  8  gives  the  corresponding  experimental  results.  At  low 
angles  of  attack,  the  computed  streamline  patterns  on  the  flat  part  of  the 
delta  wing  are  in  qualitative  agreement  with  the  experimental  data.  (Com¬ 
parsion  is  possible  only  on  the  flat  part  of  the  wing.)  At  higher  angles 
of  attack,  the  experimental  streamlines  tend  to  go  outward.  A  possible  ex¬ 
planation  of  this  effect  is  that  pressure  may  increase  on  the  flat  plate 
region  in  the  ^-direction  (pressure  would  be  highest  on  the  centerline)  . 

This  would  tend  to  make  the  streamlines  go  outward.  On  this  assumption, 
an  option  has  now  been  put  into  the  program  so  that  pressure  can  be  in¬ 
creased  proportionally  with  <f>  on  the  flat  plate  area.  In  certain  cases, 
however,  the  streamlines  will  also  develop  a  tendency  to  cross;  this  is 
attributed  to  the  use  of  an  inadequate  pressure  gradient. 

Analysis  of  NASA  pressure  data  indicates  that  at  least  in  this  case 
study,  the  pressure  did  not  increase  when  going  toward  the  centerline.  It 
is  then  realized  that  the  streamline  pattern  on  the  flat  plate  underside 
of  the  wing  depends  on  two  factors:  (1)  the  pressure  distribution  and  (2) 
the  direction  with  which  the  streamline  merges  onto  the  flat  plate  region. 

The  latter  is  particularly  affected  by  pressure  gradient  on  the  cylindrical 
region. 


12 

Bertram,  M.H.  and  P.E.  Everhart,  "An  Experimental  Study  of  the  Pressure 
and  Heat-Transfer  Distribution  on  a  70°  Sweep  Slab  Delta  Wing  in  Hypersonic 
Flow,"  National  Aeronautics  and  Space  Administration  Technical  Report  R-153 
(1963). 
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A  Newtonian  conical  flow,  flat -plate  pressure  distribution  has  been 
incorporated  into  the  program  as  an  alternative  means  of  computing  the 
pressure  distribution  in  this  area.  The  theory  here  is  taken  from  References 
13  and  14,  although  it  should  be  noted  that  this  theoretical  develop¬ 
ment  assumes  sharp  leading-edge  wings.  The  pressure  and  streamline  patterns 
computed  by  these  formulas  are  not  radically  different  from  those  previously 
discussed. 

The  correctness  of  the  heat  transfer-calculation  methods  (in  both 
the  laminar  cases  now  being  discussed  and  the  turbulent  cases  given  later) 
is  determined  by  studying  the  spherical  nose,  since  one  can  be  most  certain 
about  the  pressure  distribution  in  that  area.  The  heat  transfer  results  on 
the  spherical  nose  are  shown  in  Figure  9.  After  the  first  few  computed  time 
intervals,  both  calculation  procedures  using  the  Lees  method  give  equivalent 
results,  and  these  have  therefore  been  shown  as  a  single  line.  Results  are 
also  shown  for  the  Vaglio-Laurin  laminar- streamline  divergence  method.  The 
Vaglio-Laurin  results  using  (H  -H  )  fall  lower  than  the  Lees  results  at  the 
higher  angles  of  attack.  This  trend  is  believed  to  be  correct,  and,  in 
general,  the  agreement  is  considered  good. 

Figure  10  shows  computed  and  experimental  heat  transfer  results  for 
the  cylindrical  leading  edge  and  flat  plate  regions  at  angles  of  attack  of 
a  =  20  and  30  degrees.  Data  are  plotted  with  respect  to  a  line  normal  to 
the  leading  line,  following  the  surface  contour.  The  computed  values  are 
reasonable  and  correct  in  trend  but  slightly  higher  (more  so  for  the  20- 
degree  case)  than  the  experimental  values.  Additional  discussion  of  these 
results  will  be  presented  subsequently.  Calculations  have  been  made  at 
higher  values  of  angle  of  attack,  40  and  45  degrees;  the  basic  trend  of 
these  results  (decreasing  Stanton  number  along  a  streamline)  for  regions 
downstream  on  the  flat  plate  segments  of  the  wing  disagreed  with  experimen¬ 
tal  values  which  showed  an  approximately  constant  Stanton  number.  The  theories 


13 

Hida,  K.,  "Thickness  Effect  on  the  Force  of  Slender  Delta  Wings  in 
Hypersonic  Flow,"  American  Institute  of  Aeronautic  and  Astronautics  Journal 
New  York,  pp.  427-433  (Mar  1965). 

14 

Polak,  A.  and  T.  Y.  Li,  "Three-Dimensional  Boundary-Layer  Flow  over 
a  Flat  Delta  Wing  at  a  Moderate  Angle  of  Attack,"  American  Institute  of 
Aeronautics  and  Astronautics  Journal  New  York,  pp.  233-240  (Feb  1967). 
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presented  in  Reference  12  for  the  wing  centerline  predict  approximately  the 

same  incorrect  results  downstream.  This  discrepancy  is  believed  to  be  due 

to  differences  in  the  basic  streamline  pattern  at  these  higher  angles  of 

attack  as  already  discussed.  At  lower  angles  of  attack,  agreement  between 

computed  results  and  theory  is  satisfactory  on  the  cylindrical  region: 

however,  there  is  some  disagreement  for  the  flat  plate  part.  This  difference 

could  be  partially  explained  by  limitations  in  the  theory  itself.  As  will 

be  discussed  further,  the  theory  is  limited  to  the  nose  regions  since  the 

external  flow  has  not  yet  expanded  to  a  high  velocity  in  that  area.  (The 

2 

theory  breaks  when  u  /(H  -H  )  is  not  much  less  than  1.)  It  can  then 
'  e  e  w 

be  noted  that  the  greatest  external  flow  velocities  u^  occur  on  the  aft  part 
of  the  flat  plate  region  at  low  angles  of  attack. 

Results  presented  in  Figure  10  are  essentially  heat  transfer  values 
plotted  along  lines  normal  to  the  leading  edge.  (This  represents  a  reeval¬ 
uation  of  the  printout  data  which  are  in  the  form  of  heat  transfer  values 
along  the  streamlines.)  Plotted  in  this  manner,  the  computed  heat  transfer 
results  indicate  a  falloff  in  heat  transfer  values  when  moving  outward  along 
the  span,  constant  s/d,  on  the  cylindrical  leading  edge;  in  contrast,  the 
heat  transfer  values  seem  to  fall  approximately  on  a  single  curve  on  the 
flat  plate  region.  This  type  of  information  could  be  useful  for  design 
purposes . 

Finally,  for  the  laminar  flow  case,  calculations  were  also  performed 

0  8 

at  a  =  20  degrees,  using  the  laminar  reference  values  p  y  =(p  p  ) 

0  2  r  r  e  e  . 

(PwPw)  '  (an  option  in  the  program)  suggested  by  Vaglio-Laurin .  However, 

these  results  were  not  as  good  as  those  obtained  using  the  Eckert  reference 

values. 

For  the  case  turbulent  flow,  the  correctness  of  the  heat  transfer 
methods  was  determined  by  comparsion  (matching  data  at  the  maximum  point) 
with  hemispherical  nose  data  of  Thomas  et  al.^  All  the  various  techniques 
used  seemed  to  give  reasonable  results. 


Thomas,  A.  C.  et  al.,  "Advance  Re-Entry  Systems  Heat  Transfer  Manual 
for  Hypersonic  Flight,"  Wright-Patterson  Air  Force  Base,  Flight  Dynamics 
Laboratory  Technical  Report  65-195,  Boeing  Co.  D2-85029-1.  Contract 
AF33(657)7132  (Oct  1966). 
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The  computed  heat  transfer  results  for  the  case  of  turbulent  flow 
were  also  compared  with  NASA  experimental  values  for  the  type  of  delta  wing 
under  consideration.^  Only  a  few  sets  of  data  giving  both  pressure  dis¬ 
tribution  and  heat  transfer  are  presented,  and,  unfortunately,  the  reported 
heat  transfer  results  are  limited.  Comparsion  was  made  with  a  delta  wing 
swept  back  73  degrees  at  Mm  =  7,  a  =  20  and  30  degrees  (30  degrees  was  the 
highest  value  reported) ,  The  results  (Figure  11)  tended  to  be  somewhat 
similar  to  the  laminar  case  in  that  agreement  was  best  for  the  30-degree 
example.  As  was  found  by  Tai,^  results  are  better  for  the  modified  tech¬ 
nique  of  Rose  et  al.  than  for  the  streamwise  divergence,  turbulent  flow 
method.  Only  centerline  turbulent  flow  data  are  reported  in  Reference  16, 
whereas  the  present  program  calculates  streamlines  going  about  the  body 
(from  the  spherical  nose  to  the  cylindrical  leading  edge  to  the  flat  plate 
area)  which  eventually  intersect  the  centerline.  Because  of  this  difference 
in  calculation  procedures,  the  present  comparsion  is  actually  a  worst-case 
type  of  comparsion. 

In  the  development  of  his  theory,  Vaglio-Laurin  states  various  re¬ 
quirements  : 

1.  The  flow  must  be  hypersonic 

2.  The  wall  must  be  cold  and 

2 

3.  u  / (H  -H  )  must  be  very  much  less  than  1. 

e  ^  e  w  7 

Basically,  these  requirements  are  needed  to  ensure  that  the  boundary  layer 
flow  direction  will  be  the  same  as  the  external  flow  direction;  this  would 
not  be  the  case,  for  instance,  in  the  subsonic  or  transonic  regimes.  Re¬ 
quirement  3  indicates  that  the  validity  of  the  theory  will  decrease  as  one 
moves  downstream  of  the  nose  onto  regions  of  the  body  that  exhibit  a  low 
effective  angle  of  attack  relative  to  the  free  stream.  In  view  of  this,  the 
validity  of  the  theory  at  a  given  point  on  the  body  may  be  expected  to  in¬ 
crease  with  increasing  angle  of  attack.  The  pressure  formulas  used  in  the 
present  study,  however,  lose  validity  at  the  higher  angle  of  attack  values. 
It  is  therefore  not  unexpected  that  the  agreement  between  theory  and  data 
maximizes  at  a  moderate  value  of  angle  of  attack,  30  degrees. 


Nagel,  A.  L.  et  al.,  "Analysis  of  Hypersonic  Pressure  and  Heat  Transfer 
Tests  on  Delta  Wings  with  Laminar  and  Turbulent  Boundary  Layers,"  National 
Aeronautics  and  Space  Administration  Contractor  Report  535  Boeing  Co.  Contract 
NAS  1-4301,  (Aug  1966). 
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A  second  point  to  be  discussed  is  the  influence  of  the  shock  wave 
shape  characteristics  for  the  delta  wing.  The  important  point  here  is  that 
the  boundary  layer-edge  flow  (which  is  eventually  entrained  into  the  bound¬ 
ary  layer)  is  calculated  as  an  isentropic  expansion  of  the  flow  downstream 
of  the  shock  wave.  The  present  program  assumes  that  all  the  flow  that  is 
entrained  into  the  boundary  layer  has  passed  through  the  normal  shock  at  the 
nose;  this,  of  course,  is  an  approximation,  especially  in  the  case  of  thin 
wings  at  low  angles  of  attack.  A  second  approach  is  to  initiate  the  expansion 
across  an  oblique  shock  determined  by  the  flat  plate  part  of  the  wing.  The 
flow  is  between  the  two  extremes. 

The  feature  that  characterizes  three-dimensional  boundary  layers  is 
the  fact  that  the  streamline  path  is  to  be  evaluated  as  part  of  the  analysis. 
This  is  in  contrast  to  the  axisymmetric  case  where  the  streamline  path  is 
determined  completely  by  symmetry.  For  this  reason,  the  influence  of  pressure 
distribution  on  the  heat  transfer  characteristics  is  doubly  cumulative  in 
three-dimensional  f lows--pressure  influences  both  streamline  path  and  local 
thermodynamic  variables.  Better  pressure  distributions  should  produce  more 
accurate  heat  transfer  values. 

CONCLUSION 

A  computer  program  has  been  generated  to  calculate  the  heat  transfer 
characteristics  on  the  windward  side  of  a  delta  wing  in  hypersonic  flow. 

The  program  utilizes  the  small  crossflow  theory  of  Vaglio-Laurin  which  is 
limited  to  hypersonic  speeds,  cold  walls,  and  moderate  boundary  layer-edge 
velocities.  This  program  may  be  used  either  for  analyzing  the  thermal 
environment  on  a  proposed  configuration  or  for  studying  the  parameters  that 
can  effect  heat  transfer  for  a  three-dimensional  configuration. 

Comparsions  of  computed  results  and  experimental  delta  wing  heat 
transfer  data  appearing  in  the  literature  show  reasonable  agreement, 
especially  for  cases  where  the  limitations  of  the  theory  are  not  exceeded, 
and  the  pressure  distributions  used  as  input  are  most  accurate. 
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-ure  1  -  Cylindrical  Leading  Edge  Geo 


Figure  2  -  Lower  Wing  Surface 
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Figure  3  -  Spherical  Nose  Part  of  Delta  Wing 


P/P0 


Figure  4  -  P/PQ  on  Spherical  Nose 
(Measured  values  are  from  Reference  12) 
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Figure  6  -  Centerline  Pressure  Distribution 
(Experimental  values  are  from  Reference  12) 
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Figure  9  -  Heat  Transfer  Distribution  on  Spherical  Nose 

(Laminar  Flow) 

(Experimental  values  are  from  Reference  12) 
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Figure  10b  -  a  =  30  Degrees 
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Figure  11  -  Heat  Transfer  Distribution  on  Centerline  for  Delta  Wing 

with  Turbulent  Boundary  Layer 

(Data  from  Reference  16) 
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